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ABSTRACT 



The baryons inside high-redshift halos with virial temperatures above T ~ 10 4 K cool radiatively 
as they condense inside dark matter potential wells. We show that the release of the gravitational 
binding energy, over the halo assembly time-scale, results in a significant and detectable Lya flux. At 



the limiting line flux 



10 



-19 



erg s cm asec of the Next Generation Space Telescope, several 



sufficiently massive halos, with velocity dispersions a <; 120 km s -1 , would be visible per 4' x 4' field. 
The halos would have characteristic angular sizes of « 10", would be detectable in a broad-band survey 
out to z ~ 6 — 8, and would provide a direct probe of galaxies in the process of forming. They may be 
accompanied by He + Lya emission at the 10% level, but remain undetectable at other wavelengths. 
Our predictions are in good agreement with the recent finding of two Lya "blobs" at z — 3.1 by Steidel 
et al. (1999). 

subject headings: cosmology: theory - early universe - galaxies: evolution - galaxies: ISM 



1. INTRODUCTION 

In standard Cold Dark Matter (CDM) models for the 
formation of structure in the universe, a significant popu- 
lation of DM halos collapse at high redshift (z ^ 5). Pro- 
vided that the gas inside such a halo is heated to a virial 
temperature of T > 10 4 K (or T > 10 2 K if H 2 molecules 
are present, see, e.g., Haiman, Abel & Rees 2000), it can 
cool on a short time-scale, and continue condensing. This 
condensation process can be regarded as a cosmological 
cooling flow: a high-redshift and low-mass analogue of 
cooling flows found in clusters of galaxies (e.g. Fabian 
1994). 

The behavior of the gas in early DM condensations 
should resemble the formation of normal disk galaxies, 
which has been investigated by several authors in semi- 
analytic schemes (White & Rees 1978; Fall & Efstathiou 
1980; Mo, Mao & White 1998 [MMW] ; van den Bosch & 
Dalcanton 2000), as well as numerical simulations (Katz 
1991; Navarro & White 1994; Navarro, Frenk & White 
1997 [NFW] ; Navarro & Steinmetz 1997; Moore et al. 
1999). In all of these studies, most of the baryons in the 
halo cool and settle into a rotationally supported disk with 
an approximately exponential surface density profile. In 
the simplest picture, the disk material originates as smooth 
gas, collapsing isothcrmally from the virial radius i? V i r to 
its final orbital radius of ~ Ai? V i r , where A ~ 0.05 is the 
typical spin parameter. The numerical simulations reveal a 
more complex process, where the gas forms smaller clumps 
early on; these clumps then progressively merge together, 
collide, and dissipate, to form the larger systems. 

Irrespective of the dynamical details of the collapse, the 
disk material must dissipate a minimum energy ~ M^ishV^ i 
where Mdisk is the final mass of the disk and v c is the 
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circular velocity of the DM halo in which it is embed- 
ded. The most natural mechanism for this dissipation is 
radiation from collisionally excited H and He, as well as 
Bremsstrahlung radiation for more massive (T vlI <; 10 6 K) 
DM halos. Although this cooling radiation has been as- 
sumed to be present in studies of the formation of disk- 
galaxies, its observability has typically not been addressed. 

In this Letter, we use a simple toy-model to quantify 
the Lya fluxes expected from the baryonic material cool- 
ing and settling into the DM potential wells at high red- 
shift. We estimate the number of sources that would 
be detectable with the Next Generation Space Telescope 
(NGST), and argue that these sources could be identified 
based on their their large angular extent, relatively shal- 
low surface brightness profiles, and their low luminosity 
at other wavelengths. A detection of the cooling radiation 
would be a novel and direct way to study proto-galaxies 
in the process of their assembly. Throughout this Letter, 
we assume the background cosmology to be ACDM, with 
the parameters Q m = 0.3, Q& = 0.04, f^A = 0.7, h = 0.65, 
and erg = 1. 

2. COOLING RADIATION 

We first estimate the energy budget and total flux of the 
cooling radiation from a halo under the following simpli- 
fying assumptions: (1) initially, the baryons and the DM 
both follow the same density profile, apart from a normal- 
ization by f2;,/O m , (2) the baryons settle into a rotation- 
ally supported disk, and (3) the energy radiated emerges 
as Lya radiation. 

A parcel of gas in a massive halo with virial tempera- 
ture significantly above 10 4 K would be shock-heated to a 
high temperature, and initially cool via either He + line- 
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cooling, or via Bremsstrahlung. This initial, fast cooling 
would release the thermal energy of the gas, and reduce its 
temperature to ~ 10 4 K. 3 Since the 10 4 K gas is no longer 
pressure supported, it contracts on a dynamical timescale, 
and moves inward in the DM potential well. During this 
contraction, it is likely to remain roughly isothermal, and 
radiate mostly Lya radiation (see discussion below). 

The energy dissipated in this contracting phase must be 
approximately the gravitational binding energy. We now 
compute this energy explicitly for the gas moving in the 
NFW potential. We assume that the gas and DM are ini- 
tially both distributed according to 



P(r) = Po 



1 



{r/r s ){l + r/r s f 



(1) 



and that the final state is a rotationally supported expo- 
nential disk embedded in the DM halo. We use the solution 
for the disk surface density S(r) in MMW (their eq. 28) in 
the final configuration, which takes into account the grav- 
itational effect of the disk on the halo in a self-consistent 
fashion. A gas particle initially located at radius Ri set- 
tles into a circular orbit at a final radius R / , such that the 
initial and final enclosed gas masses are equal, 



drAirr 2 p gas (r 



dr27rr£ gas (r). 



(2) 



In the solution for E gas (r), we adopt the typical concen- 
tration index c = R v - 1T /r s = 10, spin parameter A = 0.05, 
and assume that all baryons end up in the disk. We find 
contraction factors of Ri/Rf ss 10, varying from a value 
of » 8 for particles near the the surface i? v ir to a peak of 
ps 13 near half the virial radius. Note that more realisti- 
cally, a distribution of values for A should be used, causing 
an intrinsic scatter in the luminosity of a halo of given ini- 
tial size. This could enhance the number of halos visible 
at a given Lya flux. 

The total energy released in cooling by a particle of mass 
m in this process is 



EcncA —W — -Ekin = 



(If 



R f 



GM{r)m GM(R f )m 

—y* r} ■ 



(3) 



where W is the work done on the particle by the gravita- 
tional field, -Ekin is the kinetic energy of the particle in its 
final circular orbit, and M(r) is the total mass interior to 
radius r. Summed over all particles, we have 

W = dRA^RhzM) dr^f± (4) 



Rt 



1 



dr2itrY>(r) 



GM'(r) 



(5) 



Note that in eq. [|, we have used the initial gas distribu- 
tion /9 gas and total (gas + DM) enclosed mass M(r). For 
simplicity, we have assumed that the initial NFW distri- 
bution determines both quantities. By contrast, in eq. ^ 
£(r) and M'(r) refer to the surface density of the gaseous 
disk, and the total enclosed mass (DM+disk) in the fi- 
nal equilibrium state, taken from the numerical solution 



in MMW. Under these assumptions, we find that the total 
radiated energy is 



E, 



cool 



= W - Eun = 1.9£ b 



ind i 



(6) 



where Ebind = f v " drAirr 2 p gas GM(r)/r is the total ini- 
tial gravitational binding energy. 

In an early three-dimensional smooth particle hydrody- 
namics (SPH) simulation, Katz (1991) explicitly computed 
the energy dissipated during the collapse of a 7 x 10 11 M Q 
perturbation in a standard CDM. Adopting the parame- 
ters of this halo, we find the energy radiated in our spher- 
ical approach to be 1.1 x 10 59 erg. This is only w 30% 
smaller than the value 1.5 x 10 59 obtained in the simu- 
lation, lending support for our simple estimates for the 
overall energy budget. 

The timescale over which the dissipated energy is re- 
leased is related to the assembly of the halo. In the sim- 
ple picture of monolithic collapse in a pre-existing spheri- 
cal DM halo, the initial thermal energy would be released 
quickly on a cooling time; subsequently, the gas would con- 
tract isothermally on a dynamical time. Since the typical 
overdensi ties are of order ~ 100; this would correspond 
to - 1/V^LQ0 = 10% of the Hubble time t Buh (z). How- 
ever, in the hierarchical structure formation scenarios of 
CDM cosmologies, halos are built up gradually, as a re- 
sult of collapse on smaller scales, mergers, and continuous 
accretion of fresh gas. The typical time-scale on which 
this process occurs can be estimated within the extended 
Press-Schechter "merger tree" formalism. This formalism 
yields the probability distribution of ages dp/dt(M, z, th) 
for halos of mass M existing at redshift z, where the char- 
acteristic "age" , th, is defined to be the time elapsed since 
the halo had acquired half of its present mass (Lacey & 
Cole 1993, eq. 28). 

For a halo of mass M at redshift z, we adopt the defi- 
nition 

' !p (M 7 z,t) 



dt 



dt 



0.5 



(7) 



for the characteristic time td over which the gas assem- 
bles and dissipates its energy. This definition ensures that 
half of the halos build up half of their final mass within 
the time td- The time td is a fraction of the Hubble 
time and depends on the halo mass. On average, larger, 
rarer halos are assembled more recently than typical ha- 
los corresponding to ler density perturbations. In general, 
for the halo masses relevant for our discussion, we find 
0.15 <, td/tYi u h{z) <, 0.30, slightly longer than the dynam- 
ical time. 

3. DETECTABILITY OF LYMAN a HALOS 

Armed with the energy budget and characteristic dis- 
sipation time, we can now make simple estimates for the 
expected Lya fluxes from gas cooling in halos. The total 
line flux from a halo of mass M at redshift z is simply 
given by 



F a (M,z) 



E cool (M,z)/t d (M, z) 
^df(z) 



(8) 



where di{z) is the luminosity distance to redshift z. 
A halo will subtend an angle 26 vu = 2R vlr /dp L ~ 

This could result in additional He Ly a or Bremsstrahlung cooling signals comparable to the H Lya radiation discussed here. 
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12"(i? v ; r /50 kpc)[dA(z = 5)/cIa{z)], and will be easily re- 
solved at NGSTs angular resolution of ss 0.05". 

To assess the observability of our sources, we use the 
proposed instrumental characteristics of NGST with an 
8m diameter mirror. The expected signal-to-noise ratio 
in an observation of duration t is 



S_ 

N 



AftF a (l + z)/hv a 



(AftAXnO. 
which yields approximately 



2 R 



s_ 

N 



A 



1/2 



50 m 2 
2 x 10 4 



J_ 

0.5 



1/2 



t 



h 



2o 



1/2 



M 



2.5 x 10 5 s 

4/3 



3 x 10 11 M 



1/2 



l + Z 



6 



(9) 



1/2 



(10) 



where A is the mirror area, / is a fiducial sys- 
tem throughput, R = AA/A bs is the assumed fil- 
ter width, t is the integration time, and R g = 
h. 2s(A/l. 25/im)" 15 photons cm" 2 s" 1 sr" 1 A" 1 is the 
zodiacal background, extrapolated to wavelengths shorter 
than 1.25 /Ltm from Kelsall et al. (1998). The background 
turns over at A < 0.5/im (Leinert et al. 1998), so that 
our extrapolation is an overestimate at these wavelengths 4 . 
We have approximated the angular diameter distance as 
di oc (l+z) ' 75 , and the dissipation time as id oc (l+z)" 1 - 5 . 
These scalings are accurate to within a factor of two in the 
relevant mass and redshift range; in our calculations we use 
the exact expressions. 

The expected line-widths are of order ~ 30[(1 + z)/6]A 
(see below), and a narrow band filter with a compa- 
rable width would minimize the sky background. In 
an initial blind search for the Lya halos, however, we 
find that it is more advantageous to use a broader fil- 
ter, which covers a larger redshift range. Eq. [l^ re- 
veals that the minimum detectable halo mass scales with 
the filter width approximately as M oc AA 3 / 8 ; the to- 
tal number of halos visible per observation then scales as 
N oc dzM(dN/dM) oc AA/M oc AA 5/8 . Accordingly, in 
the estimates below we adopt a broad-band filter width of 
AA/A bs = 0.2 (R — 5). The limiting surface brightness 
near the wavelength of lfim in a 2.5 x 10 5 second obser- 
vation, is F min sa 10" 19 erg s" 1 cm" 2 asec" 2 (S/N=5). In 
comparison, the limiting line fluxes in current searches for 
high— z Lya-emitters are typically « 10" erg s" 1 cm" 2 . 

In Table @, we list examples for the detectability of 
Lya halos at NGST's limiting line flux at different red- 
shifts, and summarize the properties of halos at the de- 
tection threshold. In order to produce the required flux, 
halos located at redshifts 2 < z < 8 must have circular 
velocities of around 170 km s , or total masses above 
Mmin « 10 11 Mq. Note that because of the smaller an- 
gular sizes at higher redshift, and because of the decreas- 
ing background at longer wavelengths, a smaller halo is 
detectable at higher redshifts. Interestingly, the circular 
velocity corresponding to the detection threshold is nearly 
independent of redshift. This is apparent from eq. since 
v 2 oc M 4 / 3 (l + z) 2 . In order to be detectable at z ~ 3 us- 
ing current instruments, the required circular velocity is 
larger, « 300 km s , corresponding to a total mass near 



7 x 10 12 M . In the last column of Table [jj we have shown 
the surface density of halos more massive than M m - ln , 

Nn-AzMl^* r dM^- (11) 
dzdil t Kvlh (z) JA/ min ( 2 ) dM 

where dV/dzd£l is the cosmological volume element, 
dn/dM is the mass-function of halos taken from the Press- 
Schechter (1974) formalism, td/tHub(z) is the duty-cycle, 
AVI = 4' x 4' is the proposed field of view of NGST, and 
Az = 0.2(1 + z) is the redshift range probed by our fil- 
ter. As the table shows, ~ 30 halos are detectable within 
each NGST survey field, with ~1 halo still visible from 
a redshift of z » 6 — 7. It should be kept in mind 
that much fainter halos at higher redshifts could be de- 
tectable through a narrow filter. For instance, adopting a 
width R = 1, 000 would allow a detection of a halo with 
M « 10 10 M , v c = 90 km s" 1 out to z = 10 (S/N=5, 
t = 2.5 x 10 5 s). 

At the sensitivity threshold of current instruments, the 
predicted surface density of Lya halos near z = 5 is ~ 600 
per square degree. This value is « 4% of the surface den- 
sity of Lya emitters at z = 4.5 found in a recent narrow- 
band search with the Keck Telescope (Hu et al. 1998). 
The last line of Table [l] shows the number of halos ex- 
pected at the sensitivity level of the Steidel et al. (1999) 
data. They found two Lya "blobs" , compared to our pre- 
diction of ~ 0.4 such objects in a 78 sq. amin field with 
a three-dimensional overdensity of 15, and in a redshift 
interval of Az = 0.06. Their halos have angular sizes of 
15" and 17", which are w 50% smaller than the size we 
obtained from the virial radius. 

Table 1 

Examples for the detectability and surface density of Lya halos. 



z 


A/ mi „ 


t d 


T vir 




F Q /10" iY 


e vir 


iVtot 


IO^Mq 




10 H K 


km/s 


crg/s/cm 2 


ascc 


4' X4' 


2 


7.5 


0.34 


2.0 


166 


7.96 


15.4 


11.2 


3 


5.2 


0.28 


2.1 


170 


3.97 


11.0 


9.1 


4 


3.9 


0.24 


2.2 


171 


2.41 


8.8 


5.7 


5 


3.0 


0.22 


2.2 


172 


1.63 


7.4 


3.1 


6 


2.3 


0.19 


2.2 


171 


1.16 


6.5 


1.5 


7 


1.9 


0.18 


2.1 


170 


0.87 


5.7 


0.7 


8 


1.6 


0.16 


2.1 


169 


0.69 


5.2 


0.3 


3.1 


70 


0.22 


12 


400 


3.7 


25 


0.07 



4. FURTHER CONSIDERATIONS 

4.1. Details of Collapse and Cooling 

The Lya cooling signal proposed here is independent of 
the fate of the dissipating gas in the final stages of col- 
lapse, i.e. the details of disk-formation, fragmentation, 
and star-formation. Nevertheless, a few caveats must be 
kept in mind. Three-dimensional simulations of galaxy 
formation have shown that a galactic halo is built up by 
mergers of pre-existing subclumps, rather than by contin- 
uous accretion of smooth gas. At the current resolution 
of simulations the fraction of the infalling mass in these 
clumps is w 10% (e.g. Yoshida et al. 2000). Since the 
slope of the subclump mass-function is dN/dM oc M" 1 - 8 , 



4 See the N GST exposure time calculator at http://augusta.stsci.edu for more details 
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the simulations might already have converged on the to- 
tal clumped mass. Even stronger sub-clumping, however, 
should not effect the overall energy budget: as long as 
the gas contracts and maintains its circular velocity, it has 
to dissipate the energy given in eq. |(| Furthermore, the 
gas in small sub-clumps can be stripped from the merging 
DM sub-clumps, either by the external UV background, 
or by clump-clump collisions. The baryons falling into 
the halo are then dispersed into the inter-clump medium 
and shock-heated, implying that their contraction and dis- 
sipation closely resembles the smooth accretion scenario. 
The strength of the internal shocks are determined by the 
streaming velocities, which numerical simulations find to 
be significantly smaller than the velocity dispersion of the 
halo. As a consequence, the internal shocks are not likely 
to heat the gas to temperatures above T ss 6x 10 4 K, where 
He + line cooling would start to dominate over Lya cooling. 

We have assumed that the dissipated energy is re- 
leased in Lya cooling alone. Indeed, for a pure H+He 
gas of primordial composition, this is a safe assump- 
tion, since cooling shuts off abruptly below the colli- 
sional excitation temperature of T w 10 4 K. However, 
in the presence of molecules and heavy elements, radia- 
tive cooling is possible below 10 4 K. The cooling time for 
Lya cooling at the peak of the cooling curve is given 
by tcooi = (3/2)nfc B T'vir/« 2 A coo i, where A cool sa 2 x 
10 -22 erg s _1 cm 3 and n is the baryon number density 
corresponding to an overdensity of s» 200 at the time of 
collapse. We have used T v ; r in this expression, rather than 
the gas temperature T w 1 4 K, since it is the former that 
characterizes the gravitational binding energy that is ra- 
diated away. We find that the cooling time is between 
6 — 45% of the dynamical time, td yn = l/^6nGp for the 
halos considered here. The requirement that £ coo i « idyn 
then implies that the gas temperature will be near, but 
somewhat below, the temperature T sa 10 41 K where the 
cooling function peaks. 

Now consider the cooling of metal-enriched gas. Al- 
though non-Lya cooling is possible at T < I0 4 K, the cool- 
ing function below this temperature drops sharply. As 
long as the metalicity is Z < O.IZq, molecular cooling or 
cooling by atoms heavier than He is at least a factor of 
« 1000 less efficient then cooling via the Lya peak. Since 
initially the non-Lya cooling time exceeds the dynamical 
time by a factor of 60 — 450, and t C oolAdyn oc p~ 1 ^ 2 , Lya 
cooling dominates until the density is enhanced by a fac- 
tor of 4 x 10 3 — 2 x 10 5 , or until radial contraction factors 
of 15 — 60. We therefore conclude that the bulk of the 
binding energy is likely released in Lya cooling, until the 
metalicity builds up to near-solar levels. 

Dust absorption can strongly suppress the Lya flux es- 
caping from a medium that is optically thick to Lya pho- 
tons; this is thought to be the reason why early Lya sur- 
veys did not detect proto-galaxies (see Pritchet 1994). Re- 
cently, however, Lya emitting galaxies have been found at 
high-redshift (e.g. Hu et al. 1998), as expected in models 
with lower galactic dust abundance, and inhomogeneous 
dust distribution (Haiman & Spaans 1999). Furthermore, 
the dust abundance in the early, spatially extended, col- 
lapsing phase of the high-redshift halos is likely to be sig- 
nificantly lower than inside star-forming galaxies. 

The typical Lya luminosities we derive are ~ 
10 43-44 erg s _1 . Our halos are expected to be mostly neu- 



tral, and since they are embedded in the UV background, 
they will convert a fraction 0.5 of the incident ionizing 
radiation to Lya radiation (Gould & Weinberg 1996). It 
is interesting to note that for a UV background flux of 
5 x 10~ 22 (hu/13.6 eV)- 1 ' 7 erg s" 1 cm" 2 Hz" 1 sr" 1 , the 
resulting Lya radiation would be at a level of 5 — 10% of 
our cooling signal. 

4.2. Contribution of Lya Halos to the UV Background 

In addition to the Lya flux expected from individual ha- 
los, we have computed the contribution to the extragalac- 
tic UV background from halos below the detection thresh- 
old, by summing the flux of these sources. We find the am- 
plitude for this background to be quite small, at the level 
of £ 1% of the UV background of 20±10 nW m" 2 sr" 1 
(Bernstein 1997). Note, however, that the Lya back- 
ground from collapsing halos should always dominate at 
sufficiently high redshifts, when other UV sources, such as 
stars or quasars, have not yet formed. 

4.3. Line Profiles: Shape and Interpretation 

In order to maximize the diagnostic value of observa- 
tions of the Lya halos, it is important to understand the 
frequency distribution and surface brightness profiles over 
which the Lya photons are emitted. In a naive spherical 
model, one may assume a steady emission of Lya pho- 
tons over the halo assembly time. For a purely radial 
contraction inside the logarithmic potential of an isother- 
mal sphere, the luminosity of each shell would scale ap- 
proximately as Lir) oc J r 2 dr{nk^T vlT ) /td oc J r 2 drn 3 ^ 2 oc 
logr. Most of the luminosity would arise near the outer 
shells, and the surface brightness profile would scale ap- 
proximately as J(r) oc r~ 2 out to w 30% of the virial 
radius. A similar conclusion holds for the NFW profile, 
but with J(r) flattening to I(r) oc r _1 near r = 0. 

Although speculative, these surface brightness distribu- 
tions have interesting consequences. Since they are peaked 
towards r = 0, this could potentially make the Lya halos 
more observable than our estimates in Table [l] imply. For 
example, a fraction as 1/3 of the total luminosity is re- 
leased internal to the radius i? V h_/10. Taking the total 
flux to be 1/3 of the value in eq. g, and assuming an ex- 
tended source with an angular size of v i r /lO significantly 
increases the number of detectable sources. At redshift 
z = 5, the inner ss 0.6" region of a v c = 130 km/s halo 
would still be detectable, providing w 12 candidate sources 
per NGST field. Similarly, the innermost « 0.3" region of 
a halo with v c — 120 km/s would be visible to z ~ 10. 

The depth of the DM potential well would imply a line- 
width of ~ 200 km s _1 . The line profile is, however, signif- 
icantly modified by scattering effects. Assuming that the 
halo is mostly neutral, the Lya optical depth at the center 
of the line is t = uho-qRv^ w 10 7 . For a thermal speed of 
~ 10 km s" 1 it requires a shift in frequency by an amount 
SX/X as 3 x 10~ 3 to diffuse out of the optically thick core of 
a Voigt-profile (i.e. to a frequency where r = 1). This will 
significantly broaden the line, to ~ 1, 000 km s~ 4 ; it would 
also further flatten the surface brightness profile. Depend- 
ing on the infall velocity and neutral density structure, the 
line may exhibit a characteristic "double-hump" profile. 

It would be very interesting to extract these quantities 
directly from a three-dimensional simulation in the fu- 
ture, since the angular resolution and spectral capabilities 
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of NGST should be sufficient to map the surface brightness 
distribution, and to determine the line profile, at least for 
the brightest halos. Here we simply note that the char- 
acteristic line-widths of ~ 1000 km s _1 , as well as the 
extended, and relatively shallow surface brightness pro- 
files, would be robust indicators of the initial stages of the 
collapse. 

5. DISCUSSION AND CONCLUSIONS 

We have considered the Lya flux from baryonic gas con- 
densing and cooling inside high-redshift galactic DM ha- 
los. Since the hydrogen in the universe is reionized beyond 
a redshift z ps 6, the cooling radiation would not suffer 
strong scattering in the neutral 1GM below this redshift 
(cf. Rybicki & Loeb 1999), but the halos would still appear 
extended on the sky with an angular diameters of ps 10". 
We have ignored the radiation of the initial thermal energy 
of the gas. Most of this energy is radiated in the contin- 
uum as Bremsstrahlung, but is too faint to be observable 
with CXO or XMM. However, for sufficiently metal-poor 
halos, ps 10% of the initial thermal energy could be re- 
leased by He + Lya cooling. This emission likely occurs on 
the halo-assembly time-scale, and would accompany our 
hydrogen Lya signal at a wavelength 4 times shorter, and 
with a ps 10 times lower line flux. 

It is interesting to compare the overall binding energy 
radiated in Lya photons by the pristine halos and the stel- 
lar Lya known to be produced in "normal" galaxies. The 
total specific energy radiated in Lya is E(cool) ps 2v%, 
where v c is the circular velocity. According to Kenni- 
cut (1998), 1 solar mass of gas cycled through stars re- 
leases 3.1 x 10 49 erg in Lya. Assuming that 10% of the 
gas is converted into stars, we find a specific energy of 
-E(stars) = 1.7 x 10~ 6 c 2 . As a result, for a v c — 200km/s 
halo we find i?(cool)/i?(stars) = 0.5. This relatively large 
ratio demonstrates that the total flux in the cooling radi- 
ation can be comparable to that in stellar Lya emission 



lines. 

Our predictions for the fluxes, angular sizes, characteris- 
tic line- width, surface brightness profiles, and the number 
of sources are in good agreement with the two Lya "blobs" 
discovered at z = 3.09 in a 78 sq. amin field by Steidel 
et al. (1999), which leads us to believe that examples of 
cooling Lya halos may already have been found. Our sim- 
ple spherical model, however, predicts a "double-hump" 
line profile, which is not seen by Steidel et al. Further 
work is required to assess whether the 3D structure of the 
collapsing halo can wash out the predicted double peak. 

Our results show that a few dozen cooling halos could 
be detected in the future out to z ps 6 — 8 in a single 4' x 4' 
field, observed for 2.5 ksec with NGST. In order to extend 
the redshift range over which this study is feasible, a large 
area survey with NGST would be desirable, with a rela- 
tively broad filter (to capture rarer, higher- z halos, over a 
larger redshift range), as well as wavelength coverage ex- 
tending down to ps 0.4^m (to capture halos at redshifts as 
low as z ps 3). Using a narrow-band filter that minimizes 
the background would allow the detection of additional, 
fainter sources, and a detailed study of the brighter ha- 
los. The cooling halos would be a novel and direct probe 
of galaxy formation in the high-redshift universe, and of 
their evolutionary history from z ps 6 — 8 to lower redshifts. 
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